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Abstract 

We study a wide class of finite-dimensional su(/u|/i)-supersymmetric models closely related to 
the representations of the Yangian Y(g\(m\n)) labeled by border strips. We quantitatively analyze 
the degree of degeneracy of these models arising from their Yangian invariance, measured by 
the average degeneracy of the spectrum. We compute in closed form the minimum average 
degeneracy of any such model, and show that in the non-supersymmetric case it can be expressed 
in terms of generalized Fibonacci numbers. Using several properties of these numbers, we show 
that (except in the simpler su( 111) case) the minimum average degeneracy grows exponentially 
with the number of spins. We apply our results to several well-known spin chains of Haldane- 
Shastry type, quantitatively showing that their degree of degeneracy is much higher than expected 
for a generic Yangian-invariant spin model. Finally, we show that the set of distinct levels of a 
Yangian-invariant spin model is described by an effective model of quasi-particles. We study 
this effective model, discussing its connections to one-dimensional anyons and properties of 
generalized Fibonacci numbers. 
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1. Introduction 

In this paper we shall consider a general class of finite-dimensional quantum models which 
by construction are invariant under the Yangian algebra Y(g\(m\n)). The key to the Yangian 
invariance of these models is their connection to certain finite-dimensional representations of 
T(gl(m|»)) labeled by a type of skew Young diagrams, the so called border strips (DEI. These 
representations play a fundamental role in the study of the integrable two-dimensional conformal 
field theory related to the latter models, namely the su{m\n) WZNW model at level 1 EH3. 
The models we shall be interested in, which we shall refer to as Yangian-invariant su(;n|n) spin 
models, are characterized by the fact that their Hilbert space is a direct sum of the irreducible 
representations of Y(g\(m\n)) labeled by border strips with exactly N boxes, where N is the 
number of spins. The simplest examples of these models are certain integrable spin chains with 
long range interactions invariant under the Yangian for a finite number of spins, namely the 
su(m|n) supersymmetric versions of the Haldane-Shastry (6j-[8) and Polychronakos-Frahm |9J- 
im chains; see, e.g., (UEHsnnma. 
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As is well known, the Haldane-Shastry (HS) spin chain is a circular (translationally invariant) 
chain with equally spaced sites, the interactions between any two spins being inversely propor¬ 
tional to the square of their chord distance. This model appears in connection with a wide variety 
of topics in theoretical and mathematical physics, including one-dimensional anyons GUSHSHHl, 
conformal field theory If7374l8l , quantum chaos if79)4221 . quantum information theory |23| , and 
quantum integrability El GIB A distinctive feature of the HS chain is the fact that it can be 
obtained from a dynamical model, namely the spin Sutherland (trigonometric) model 12674281 . 
in the strong coupling limit 0. By applying the same procedure to the spin Calogero (ratio¬ 
nal) l29l [30 j and Inozemtsev (hyperbolic) m models, one obtains the Polychronakos-Frahm 
(PF) and Frahm-Inozemtsev (FI) Il32l spin chains, which share many properties with the origi¬ 
nal HS chain. In fact, although the Yangian symmetry of the FI chain has not been rigorously 
established, we shall see in Section[2]that it is isospectral to a Yangian-invariant spin model. 

It is clear by its very definition that the class of Yangian-invariant spin models encompasses 
a wide range of quantum systems. For instance, apart from the integrable spin chains mentioned 
above, it includes the family of one-dimensional vertex models studied in Refs. El ED. A 
common feature shared by all Yangian-invariant spin models is the high degeneracy of their 
spectrum, stemming from the Yangian symmetry. This statement, however, is far from precise, 
and does not shed any light on whether a certain model in this class has additional degeneracy 
due to its particular form. The main aim of this paper is precisely to perform a quantitative 
analysis of the degeneracy inherent to all Yangian-invariant spin models. To this end, we shall 
use as a concrete measure of the degree of degeneracy of a finite-dimensional quantum system its 
average degeneracy, defined as the quotient of the dimension of its Hilbert space by the number 
of distinct energy levels. We shall derive closed-form expressions for the average degeneracy 
of a “generic” Yangian-invariant su(m\n) spin model, both in the supersymmetric (mn + 0) and 
non-supersymmetric cases. In fact, these expressions provide a lower bound for the average 
degeneracy of any Yangian-invariant spin model, and thus can be used as a practical test for 
ruling out that a particular quantum system belongs to this class. 

As it turns out, the behavior of this lower bound (which we have termed “minimum aver¬ 
age degeneracy”) is rather different in the supersymmetric and non-supersymmetric cases. In 
the former case, the minimum average degeneracy is given by a simple power law in terms of 
the dimension m + n of the one-particle Hilbert space (cf. Eq. Far more interestingly, in 

the non-supersymmetric case the minimum average degeneracy is expressed in terms of gener¬ 
alized Fibonacci numbers |34| . which reduce to the standard ones for spin 1/2 (in or n equal to 
2). Using standard properties of the generalized Fibonacci numbers, we derive the asymptotic 
behavior as the number of spins tends to infinity of the minimum average degeneracy in the 
non-supersymmetric case. We find that the leading behavior is again a power law involving the 
largest root (in module) of the characteristic polynomial of the generalized Fibonacci numbers. 

Another goal of this work is to ascertain to what degree the average degeneracy of the three 
spin chains of Haldane-Shastry type (i.e., the HS, PF and FI chains) behaves as expected for 
a generic Yangian-invariant spin model. Contrary to popular belief, we find out that this is 
actually not the case. More precisely, we show that the number of distinct levels of these chains 
grows polynomially with the number of spins, whereas for generic Yangian-invariant spin models 
this growth is exponential. As explained in detail in Section [5] the ultimate reason for this 
different behavior is the fact that spin chains of HS type possess a polynomial dispersion relation. 
As a matter of fact, our results apply as well to the whole family of vertex models studied in 
Refs. 141311331. which also have this property. Moreover, when the coefficients of a polynomial 
dispersion relation are rational numbers (as is the case for the HS and PF chains, and for the FI 
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chain when its parameter is rational) we show that the number of distinct levels behaves as N k+] , 
where k is the total degree of the dispersion relation. 

Another aspect of Yangian-invariant spin models that we shall analyze in this paper is their 
connection with one-dimensional anyons. In fact, it is well known that spin chains of HS type 
provide one of the simplest realizations of anyons in one dimension via Haldane’s fractional 
statistics d. For instance, as shown in the latter reference, the spectrum of the su(2) HS chain 
is the same as that of a system of spin 1/2 anyons with exclusion parameter g - 1/2. This result 
was essentially extended to the su(m) case in (5), with g = 1/m. In this paper we focus instead 
on the set of distinct energy levels of an su(/«) Yangian-invariant spin model, and show that it can 
be obtained from the spectrum of an effective model of quasi-particles corresponding to the l’s 
in a Haldane motif 0. In particular, in the su(2) case this effective model is simply a system of 
anyons with g = 2. In the general su(m) case we compute in closed form the statistical weights 
of these effective models, and use them to derive in a novel way several combinatorial identities 
for generalized Fibonacci numbers. 

This paper is organized as follows. In Section [2] we review the description of the repre¬ 
sentations of the Yangian Y(g\(m\n)) labeled by border strips, and give a precise definition of 
Yangian-invariant su(m|n) spin models. We also recall the definitions of the su(m|n) spin chains 
of HS type, and discuss their spectra and dispersion relations. Section[3]is devoted to the com¬ 
putation of the minimum average degeneracy of an arbitrary Yangian-invariant spin model, in 
both the supersymmetric and non-supersymmetric cases. We apply our results to the su( m\n) 
supersymmetric version of Inozemtsev’s elliptic chain l35l . showing that in the su(2), su(3) and 
su(2|l) (or equivalently su(l|2)) cases it is not a Yangian invariant spin model for a wide range of 
values of N. Using several properties of the generalized Fibonacci numbers, we also determine 
the asymptotic behavior of the minimum average degeneracy of a Yangian-invariant spin model. 
In Section [4] we study the average degeneracy of Yangian-invariant spin models that are also 
translationally invariant, like the Haldane-Shastry chain. As an example, we consider the su( 111) 
Inozemtsev chain, which is known to be (isospectral to) a Yangian- and translationally-invariant 
spin model l36l . As mentioned above, in Section[5]we present a detailed analysis of the average 
degeneracy of spin chains of HS type and, more generally, of Yangian-invariant spin models with 
a polynomial dispersion relation. In Section [6] we introduce the effective models describing the 
set of distinct levels of a non-supersymmetric Yangian-invariant spin model, and discuss their 
interpretation in terms of anyons and their connection with several identities satisfied by gener¬ 
alized Fibonacci numbers. The paper ends with a brief section in which we summarize our main 
results. 

2. Yangian-invariant spin models 

In this section we shall provide a precise definition of the class of Yangian-invariant models 
on which we shall focus in this paper. For the reader’s convenience, we shall begin by recalling 
the precise definition of the Yangian Y(g\(m\n)). Following the original paper |[37l and Ref. f38l . 
we first give an explicit matrix realization of the Lie superalgebra gl(m\n). To this end, let us 
introduce a Z 2 grading p in the vector space C" !+ " by setting p{v) = 0 if v e C m x {0) and 
p(v) — 1 if v e {0) X C". This induces a natural grading in the space M m+n ( C) of (m + n) x (m + n ) 
complex matrices given by 

deg(F^) = p(a) + p(/J ), 

where p{a) = p(e a ), {e\,... , e m+n ) being the canonical basis of C m+ ", and E a P denotes the ma¬ 
trix whose only nonzero entry is a 1 in the a-th row and /3-th column. By definition, the Lie 
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superalgebra gl(m\n) is the vector space M,„ W ,(C) with the the graded commutator defined by the 
customary formula 

[A, B] ± = AB - (_l)degA+degB BA 

on matrices with well-defined degree, and extended by linearity to the whole space. The Yan- 
gian Y(g\{m\n)) is then defined by means of the solution of the Yang-Baxter equation given by 
the /^-matrix 

R(u ) = u - hP , (1) 

where u is the spectral parameter, h e C and P is the graded permutation operator defined by 

P(e a ® ep) = (~l) p(a)m ep ® e a . 

More precisely, Y(g\(m\n)) is the associative (complex) superalgebra with generators 
T <k p , a,p = \ ,m + n, k= 1,2,..., 

defined as follows. Let 

m+n 

T(u) = T a p(u)E a/i e M m+n (Y(g\(m\n ))), 

a,f3= 1 

where the matrix elements T a p{u) are defined by 


T a p(u) + ■ (2) 

k= 1 

We also introduce a Z 2 grading in F(gl(m|n)) by setting 

de s(0 = p(“) + » 

so that the matrix 7T«) and the identity matrix II = E aa are both even. Given two even 
(m + n) x (m + n) matrices A, B with matrix elements in F(gl(m|n)), we define their graded tensor 
product by 

(A ® B)ip kl = (-\f Pi+Pk)Pi A ik B kl . 

The defining relation for the generators T® are then given by the equation 
R(u - v) l T{u) 2 T(v) = 2 T(v) l T(u)R(u - v), 

where the matrix R is defined by Eq. Q and we have used the standard notation 1 T{u) = T(ii)® IT, 
2 T(u) = II ® T(it). Using the above definition of graded tensor product, it is easy to check 
that the latter equation is equivalent to the following system for the currents T a p(u) of the Yan- 
gian Y(g\(m\n)y. 

(u - v)[T a p(u), T 7 s(v)] ± = (-\) p ‘' p ^ p P p " p ? h ■ (r y p{u)T aS (v) - T y p(y)T aS (u )) ; 
a ,p ,~y ,6 = 1 ,... ,m + n . 

This system obviously determines the defining relations of the generators of the Yangian 
superalgebra Y(g\(m\n)) through Eq. ([2]). In fact, it is well known that this algebra is actually 
generated by its lowest two generators T^j and T (2 j. 
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As mentioned in the Introduction, the models we shall deal with in this paper are associated 
to finite-dimensional representations of the Yangian Y(g\(m\n)) labeled by a particular class of 
skew Young diagrams, namely border strips, that we shall describe next. We shall start with a 
few preliminary definitions, following for the most part the notation in Ref. JT] . A skew Young 
diagram A/pis the set of boxes obtained by subtracting a Young diagram p from a larger Young 
diagram A D p. Such a diagram A/p is connected if for any pair of boxes a,b e A/p there is a 
sequence of adjacent boxes starting with a and ending with h. A border strip is a connected skew 
Young diagram not containing any 2x2 block of boxes. In other words, a border strip is a collec¬ 
tion of boxes arranged in columns, such that the top box in each column is adjacent to the bottom 
box in the column to its right. Thus a border strip is determined by a vector k = (ki ,..., k r ), 
where k, is the number of boxes in the ;-th column, counting from right to left; cf. Fig. [T] For 
this reason, from now on we shall usually identify a border strip with its corresponding vector of 
column lengths k. An su(m\n) semi-standard tableau (SST) is a particular filling of the boxes in 



ki 


t 



t 


Figure 1: Border strip (k\ ,k 2 ,..., k r ). 

a border strip k with integers in the range —n, -n + 1 1 according to the following two 

rules: 

1. The integers in each column (read from top to bottom) and each row (read from left to 
right) form a nondecreasing sequence. 

2. The nonnegative (respectively negative) integers in each row (resp. column) form an in¬ 
creasing sequence. 

Remark 1. The use of the term su(m|«) in the previous definition is motivated by the fact that, 
as we shall explain in detail below, these tableaux appear in a natural way in the description of 
the spectrum of several well known spin chains possessing su(m|«) supersymmetry like, e.g., the 
Haldane-Shastry chain. 

We can obviously identify an su(m|n) SST with the sequence 

s = (si. sn) , Si e {-n, —n + 1,..., m - 1}, 

of its integers read from top to bottom and from right to left; we shall call such a sequence the 
tableau’s spin configuration. For instance, the tableau in Fig.[2]is an su(2|3) (or more generally 
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su(/n|«), with m ^ 2 and n ^ 3) SST with spin configuration s = (-3,1,1,0,-2,-1,-1), 
associated to the border strip k = (3,1,2,1). We shall write s e k to indicate that s is the spin 
configuration of an su(m|n) SST whose underlying border strip is k. We shall also say that two 
spin configurations s, s' are equivalent, and use the notation s ~ s', provided that s, s' ek for a 
certain border strip k. 





-3 




1 


-2 

0 

1 

-1 

-1 



Figure 2: su(2|3) SST with spin configuration s = (-3,1,1,0, -2, -1,-1). 

As is well-known, the set of all border strips is in one-to-one correspondence with an im¬ 
portant class of finite-dimensional irreducible representations of the Yangian Y(g\(m\n)) whose 
carrier spaces (the so called tame modules ) are irreducible under the action of a distinguished 
maximal abelian subalgebra of Y(g\(m\n)) fT] [2]. Given a border strip k, the dimension of the 
corresponding irreducible representation is the number of all possible su(n?|n) SST associated 
with k. It is clear that any vector s e {—n,—n + 1 ,...,m— 1}^ can be regarded as the spin con¬ 
figuration of the su(n?|n) tableau constructed by placing the first component si at the top of the 
first (rightmost) column, and then placing each component s, in the box below (respectively to 
the left) of the (i - l)-th box if s, > s, \, or ,v, = Sj \ ^ 0 (resp. s, < Sj_ \, or s,- = s^i < 0). It 
immediately follows that the number of all su(»i|/i) SST with N boxes coincides with the number 
(m + n) N of all spin configurations s e {—n, -n + l,... ,m - 1} N . This implies that 

dim Vk(m\n) = (m + n ) N , 

kCP N 

where Vk(m\n) denotes an irreducible T(gl(;;j|n))-module associated with the border strip k and 
P N is the set of all partitions of the integer N (with order taken into account). Thus the direct sum 
©kefv ^k( m \n) has the same dimension as the Hilbert space of a spin chain with N sites and two 
species of particles, one with m and the other with n internal degrees of freedom. In fact, there 
is an important class of Yangian-invariant models whose Hilbert space decomposes as the direct 
sum (J) ke?> Vk(m\n), namely the su(m |n)- supersymmetric Haldane-Shastry and Polychronakos- 
Frahm chains IfTl l4l [TTl fT2l 1201 l39l that we shall now briefly describe. 

To begin with, the Hilbert space of these chains is spanned by the canonical spin basis 

|si) ® • ® |sjv) = |si,..., sjv) = |s), Si 6 {-n, —n + 1,... ,m - 1}, 

where the m bosonic (resp. n fermionic) internal degrees of freedom of the ;-th spin correspond 
to nonnegative (resp. negative) values of Sj. The action of the su(m|/;)-supersymmetric spin 
permutation operators S (with i < j) on the canonical spin basis is then defined by 

s‘'” w \si,...,Si,...,Sj,..., S N > = (T ij(s) \Si,...,Sj,...,Si,...,S N ), 

with <Tjj( s) = -1 when s, and sj are both fermionic, or when ,y, and sj are of opposite type and the 
number of fermionic spins occupying the sites i + 1,..., j - 1 is odd. In terms of these operators, 
the Hamiltonians of the su(»;|n) HS and PF chains can be taken as |f20l[40l 

h = X 

1 ^i<KN 
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( 3 ) 












where 



\ sin 2 (tt(/ - j)/N ), 


for the HS chain, 
for the PF chain. 


(4) 


and £1 < • ■ ■ < £,v are the zeros of the Hermite polynomial of degree N. 

For the HS and PF chains, there is a simple explicit formula for the energy associated with the 
irreducible representation labeled by a border strip k e P,y that we shall now recall. To this end, 
it is convenient to introduce an alternative notation for the border strips in terms of the so-called 
motifs 0119 ■ More precisely, to a given a border strip k e P,y we shall associate the vector 
8 e { 0 , l }^" 1 with components 



1, i - k\, k\ + fc 2 , ... , k\ 4 -+ k r — i, 

0 , otherwise. 


We shall call this vector 8 the representing the border strip k, and shall refer to the posi 


tions Kj = k\ + -H h of its nonzero components as its rapidities. For instance, for the border 

strip k = (3,1,2,1) e P 7 in Fig. [2] the rapidities are 3,4,6, and thus the corresponding motif is 
(0,0,1,1,0,1). For a given su(m\n) SST with spin configuration s, we define 8( s) as the motif 8 
of the border strip determined by s. It is immediate to convince oneself that 



1 , if s i+ i < St or si = s i+ i < 0 , 
0 , otherwise. 


With these definitions, the spectrum of the su(m|n) HS and PF chains with N sites (with the 
correct degeneracy for each level) can be shown to be the set of numbers nuns 


JV-l 



(5) 


with dispersion relation 



j(N - j ), for the HS chain. 


( 6 ) 


for the PF chain. 


In fact, a similar result holds for the (purely bosonic or fermionic) su(/h) Frahm-Inozemtsev (FI) 
chain (32), whose Hamiltonian is of the form Q with couplings 



Here e 2 ^ denotes the k-th root of the generalized Laguerre polynomial L" 1 , and a > 0 is a free 
parameter. In this case the spectrum is also given by Eq. (|5), with dispersion relation PD 


enU) = j(u + j~ 1 ) • 


(7) 


1 In the original definition of Haldane, the motif is the sequence of the components of S with a zero added at both 

ends. 
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It can be shown that this formula for the spectrum holds as well in the su(m\n) case. In view of 
these facts, it is natural to conjecture that the su(m|n) FI chain is also invariant under the Yangian 
Y(g\(m\n)), although to the best of our knowledge this result has not been rigorously proved. 

We shall study in this paper the degeneracy of the spectrum of the class of Y(g\(m\n))- 
invariant quantum systems whose Hilbert space 7Y V decomposes as the direct sum 

*H n = ® V k (m\n) (8) 

keP N 

of the irreducible representations of Y(g\(m\n)) labeled by border strips with exactly N boxes. 
We shall call any such system a Yangian-invariant su(m|«) spin model with N sites. As we have 
just seen, this class includes the su(w|n)-supersymmetric HS and PF long-range spin chains. The 
Yangian invariance of these spin models and the decomposition ([ 8 ]) imply that their spectrum is 
of the form 

£>n{ s), s e {-n, -n + 1 ,...,m — 1} W , 
where & N is a real-valued function satisfying 

S ~ s' => & N ( s) = & N (s'). 

Indeed, two spin configurations are equivalent if and only if they determine the same border strip 
k, which labels a unique irreducible representation V k (m\n ) of the Yangian. Since obviously 

s ~ s' <=> 6( s) = 5(s'), 

we must have 

<5,v(s) = E n (6( s)) . 

Thus the spectrum of a Yangian-invariant spin model can be equivalently described by the func¬ 
tion E n : {0, l }^ -1 —> R as the set of numbers 

E n (6( s)), S€{-n,-n + l,...,m-l} N . (9) 

We shall call E N the energy function of the system. It follows from Eq. <[9]» that the Hamiltonian 
of a Yangian-invariant spin model can be expressed as 

H=Yj e n(S)J] isXsl, 

k eP N sek 

where {|s) | s 6 k} is any orthonormal basis of the subspace V k im\n) and 6 denotes the motif 
corresponding to the border strip k. For instance, for the su(w|«) HS and PF chains the energy 
function is the linear functional 

v-1 

E N (S) = Y j e N U)S j . (10) 

j=i 

with fe'v given by Eq. | 6 |. In fact, more general Yangian-invariant spin models having a linear 
energy function ( | l()[ i with polynomial dispersion relation ev(./j have recently been studied in 
Refs. 0133). 

We shall take Eq. (J9]i as the basis of our analysis of the spectrum of a Yangian-invariant spin 
model. For this reason, our results will also apply to models like the FI chain, whose Yangian 
invariance has not been established but which is known to possess a spectrum of the form 0 
(with En given by 0 -®). 
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3. Average degeneracy and generalized Fibonacci numbers 


3.1. Minimum average degeneracy 

The main goal of this section is to compute a lower bound for the average degeneracy of the 
spectrum 0 of a Yangian-invariant spin model with N sites, defined by 


flfjV - 


(m + n) N 
f'N 


(ID 


where is the number of distinct energies. This is of course equivalent to finding an upper 
bound on Cm, a problem which we shall now address. The key observation in this respect is to 
note that, by Eq. Cm is obviously equal to the number of distinct values taken by the energy 
function Em(8), where 6 ranges over the set 


A N (m\n) = j£(s) | s e {-«, -n + l,... ,m - if]. (12) 


of all valid su(m\n) motifs with N - 1 components. Thus 


Cn < v N (m\n ), 


(13) 


where v N (m\n) denotes the cardinal of the set A N (m\n), and therefore 


<7/V ^ <7\'. m i n i 

where the minimum average degeneracy dM. min is given by 

(in + ri) N 




v N (m\n) 


(14) 


(15) 


In fact, the minimum average degeneracy is strictly less than the average degeneracy if and only 
if the energy function Em is not injective; we shall say in this case that the Yangian-invariant spin 
model exhibits accidental degeneracy. We shall also say that a Yangian-invariant spin model is 
generic if it has no accidental degeneracy. 

As we have just seen, in order to evaluate dM, m i n in closed form we only need to compute 
the cardinal of the set ( p~2] >. To this end, we shall next recall a simple characterization of all 
allowed su(m\n) motifs which shall be of fundamental importance in what follows. Consider, to 
begin with, the genuinely supersymmetric case mn 4 0. It is easy to realize that in this case any 
sequence 6 e {0, 1 } N ~ 1 is a valid motif, so that 


vat (/ m |«) = 2 n 1 , mn 4 0 . 


(16) 


Indeed, if 8 e {0,1 ) iV 1 let K\, ..., K,_ ] be the positions of the nonzero components of 8, and 
define k, = Kj - K, \ (with K {) = 0 and K r = N). Then the vector 


s = (0,...,0,-l,0,...,0,...,-l,0,...,0) 

k\ k 2 k r 


is the spin configuration of an su(w|n) SST whose corresponding motif is 8. Thus in this case we 
just have 
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mn 4 0 . 


(17) 







Consider next the purely bosonic su(m|0) case. It is clear that an su(m|0) motif cannot con¬ 
tain a sequence of m or more consecutive 1 ’s, since the spin configuration of any tableau asso¬ 
ciated to such a motif would contain a sequence of m + 1 or more distinct integers in the range 
0,1 1. Conversely, suppose that the vector 8 e {0,1 } V1 contains no sequence of m or 

more consecutive l’s. Such a 8 can be written as a succession of sequences of consecutive 0’s 
and 1 ’s of the form 

<s = (o,...,o,i,...,i,o,...) 

h Hi 

with ni ^ m - 1. It is immediate to check that the spin configuration 

s = (m - 1 ,...,m — 1 ,m — 2 ,...,m — n\,m — n\ — 1 ,m — 1 ,...), 

Zi+l 

where each entry is in the range 0,1 ,..., m — 1 on account of the condition «, ^ m-1 , determines 
an su( m) SST whose corresponding motif is 8. This shows that A ; v(/n|0) is the set of all vectors 
8 e {0, l}^ 1 containing no sequences of m or more consecutive 1 ’s. In exactly the same way it is 
proved that in the purely fermionic case Aai(0|u) is the set of all vectors 8 e {0, l}^ -1 containing 
no sequences of n or more consecutive 0’s. It immediately follows from these two facts that the 
sets A/v(m|0) and A ; y(0|mj are mapped bijectively into each other by the duality transformation 

8 — (Si,, 6 n- i) t Ay(///|0) i—> S =(1— (5i,...,l — Sjv-i') € A/v(0|m). (18) 

Hence 

v N (m\0) = v N (0\m), (19) 

so that from now on we shall restrict ourselves without loss of generality to the purely bosonic 
case. 

Remark 2. The duality transformation ( p~8| > actually arises from a mapping defined on spin con¬ 
figurations Il38ll . valid also in the general supersymmetric case. More precisely, given an su(m|n) 
spin configuration s define the vector s' = (.s',..., s' N ) by 

s'j = -s, - 1, i = 1,... ,N. 

It is immediate to check that s' is a valid su(n|w) spin configuration, and that 

tf(s') = S( s)', (20) 

where S’ is defined by Eq. m- Consider two su(m\n) and su(n|w) spin models with the same 
energy function E^. It is straightforward to show that the spectra of both models are mapped in 
a one-to-one way by the duality transformation 

E n (6) ^ E n (6') . (21) 

Indeed, it is obvious from Eq. for the spectrum that the energies of both models are in a 
one-to-one correspondence under the transformation ( |2T| . Their respective degeneracies also 
coincide, since for a given su(;n|n) motif 6 we have 


{s e {-«,-« + 1m — 1) | 8( s) = tf) = {s' e + 1,... ,n - 1} | S(s') = 8’}, 

10 



on account of Eq. ( |20| ). In particular, for the linear energy function ( |TT)| the relation between the 
corresponding spectra of two su(m\n) and su(n\m) models is simply given by 

E^ m \d’) =E 0 - E^ n \S ), (22) 

where 

N -1 

Eo = J]£N(j) (23) 

j= i 

is the maximum energy of the model(s) with a nonvanishing number of fermionic degrees of 
freedom. 

In view of the above discussion, we need only compute the cardinal v^(m) = v^(m\Q) of the 
set of su(m) = su(m|0) (purely bosonic) motifs. The key observation in this respect is that an 
su(//z) motif 6 of length N — 1 ^ m must necessarily be of the form 

k — 1 ,... ,m, 

. This observation immediately leads to the 

N ^ m + 1. (24a) 

k= l 

On the other hand, for su {in) motifs of length N - 1 < m the restriction on the number of 
consecutive l’s is vacuous, so that 

A, v (m|0) = {0, 1 }' v_l , N=l,...,m, 

and therefore 

v N (m) = 2 n ~' , N = 1,..., m . (24b) 

The recursion relation ( |24a| ), together with the initial conditions ( |24h[ >, uniquely determines 
v N (m) for arbitrary N. We shall next show that v N (m) can in fact be expressed in terms of the 
so-called m-generalized Fibonacci (in short, m-nacci ) numbers F„, defined |34| as the unique 
solution of the recursion relation 


( tf w , 0,1 . 1 ), 

k -1 

where 6 (k> is a valid su(m) motif of length N — k — 
recursion relation 

m 

v N (m) = ^ v N - k (m), 


F { n m) = Yj F TX’ ( 25a ) 

k= 1 

with the initial conditions 

F { f = ■■■ = F (m \ = 0, F (m \ = 1. (25b) 

0 m -2 ’ m -1 v ' 

Note, in particular, that for m = 2 Eqs. ([25} yield the standard Fibonacci numbers. Comparing 
Eqs. ( |24} with ([25} it is straightforward to prove that 

y N (m) = F^l m _ x . (26) 
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Indeed, it is obvious that the shifted /ra-nacci sequence F^ m _ x satisfies the recursion relation ( |24a[ >. 
Hence it suffices to show that ( |26} holds for N = 1,... ,m. This is obviously true for N = 1. 
Assuming that ( [26] ) also holds for N = 1,..., k ^ m — 1, by Eqs. ( |25} we have 

F( :l = F [ m> +■ ■ • +F \:i i+2 F ^j = 1+ £ v j+ i('«) = 1 + Z 2i = 2 * = v * +i(m) - 

7=0 7=0 7=0 

so that ( |26} also holds for N = k + 1. By induction, the latter equation is true for N = 1,..., m, 
thus establishing our claim. From Eq. GD with n = 0 we conclude that the minimum average 
degeneracy of a Yangian-invariant bosonic su (m) spin model is given by 

m N 

d N ,mm = -j-j ■ ( 27 ) 

N+m -1 

By Eq. ([T9]), the same is true for a purely fermionic su (m) spin model. 

The inequality (|TA} provides a very simple necessary condition for a quantum system to be a 
Yangian-invariant spin model, as we shall see in the following example. 


Example 1. Consider the su(m|«) supersymmetric version of the Inozemtsev elliptic chain l35l . 
with Hamiltonian 




1 

2 


Z 


l _ S M n) 
sn 2 (2 (i-j)§) 


(28) 


In the latter equation sn denotes Jacobi’s elliptic sine with real period 4 K, where K = K(k ) is the 
complete elliptic integral of the first kind 


r /2 dx 

K(k) = (29) 

V 1 - k 2 sin 2 x 

and k e [0,1) is the elliptic modulus. As is well-known, for k = 0 the latter model reduces 
to the su(m|n) HS chain (]3])-(|4| which, as we have seen in the previous section, is a Yangian- 
invariant spin model. On the other hand, for k —> 1 the Hamiltonian ( |28| ) is related to the su(m\n) 
Heisenberg chain, which does not possess Yangian invariance for a finite-number of sites. For 
k e (0,1), it is generally believed that the Hamiltonian ( [28} is not invariant under K(gl(;;;|/j)j for 
any finite value of N. We have numerically computed the average degeneracy of the spectrum 
of the model ([28} with su(2), su(3) and su(2|l) spir0(for 6 f N f 18 in the first case and 
6 f N f 12 in the remaining ones). It turns our that for k e (0,1) the average degeneracy is 
practically independent of k, so that we shall restrict ourselves to the case k 2 = 1/2. As can 
be seen from Fig. [3] for N f X the average degeneracy of the model ([28} in the cases studied 
is clearly lower than the minimum average degeneracy of a Yangian-invariant spin model (cf. 
Eqs. or ([27}). This shows that in these cases the chain ([28} is not a Yangian-invariant spin 
model (see also Ref. 11421 for the su(2) case). 


2 As shown in Ref. [38], the operators and -S^ m) are unitarily equivalent, so that the spectra of and 

ditf er by an additive constant. In particular, the average degeneracies of the su(m\n) and su(«|m) models < |28] 
coincide. 
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Figure 3: Left: logarithmic plot of the average degeneracy of the su(m) elliptic chain with k 2 = 1/2 and m — 2 (blue 
circles) and m = 3 (blue squares), compared to the minimum average degeneracy of an su(m) Yangian-invariant spin 
model with m = 2 (red circles) and m = 3 (red squares). Right: analogous plot for the su(2| 1) elliptic chain with k 2 = 1/2 
(blue circles) and the su(2|l) minimum average degeneracy (red circles). 


3.2. Asymptotic behavior of the minimum average degeneracy 

We shall next investigate the asymptotic behavior as N tends to infinity of the minimum 
average degeneracy ( [27] ) of an su(ra) Yangian-invariant spin model. To this end, note that ( |25a| ) 
is a linear recursion relation with constant coefficients and characteristic polynomial 

m 

Pm u)^r-Yj' r ~ k ■ < 3 °) 

k= 1 

It is shown in Ref. E4l that all the roots A { ' n) (i = 1,..., m) of p m (t) are simple, and that moreover 
they can be labeled in such a way that 

|d| m) | < 1, i = 1,... ,m - 1; 1 < 4" 0 = A >n < 2. (31) 

From the simple character of the roots of the characteristic polynomial it follows that the general 
solution of the recursion relation ( |25a| ) is of the form 

m 

c[ m) , « = 0,1,... , (32) 

1=1 

where the coefficients c" n are complex constants (independent of N). In order to obtain the m- 
nacci numbers, the coefficients c ' n> must be chosen in such a way that ( [32] ) satisfies the initial 
conditions ( |25b| >. i.e., 

m 

2 7 ’( 7 ') =s n,m- 1, n = 0,... ,m - 1 . 

1=1 

Note that the coefficient matrix of this linear system is the Vandermonde matrix determined by 
the m distinct roots of the characteristic polynomial p m (A), and therefore the system has always 
a unique solution. In fact, it is straightforward to show that 

c ( f = Y] (^l m) - 7T 1 . 

j;j*‘ 
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and therefore 


F (m) 
1 n 


zwTnw' 


>'= i j;j*‘ 



From this equation and Eq. © it follows that as n —> oo we have 


F™ = c m A n m + OCe-*""), 


(33) 


(34) 


with 

m —1 

c m -cw=n(^-4 m) r 

;=i 

and 

/<■„, = - min log |4 m, | > 0. 

Thus as N —> oo the minimum average degeneracy behaves as 

d N ,mm = 0 + > 

where the coefficient y m (independent of N) is given by 


Tm 


1 


C w /l 


m— 1 
m 


m -1 


n> 


A (m) ] 

i 

dm 


(35) 


(36) 


Since A m < 2 (cf. Eq. (|3T)), it follows from Eq. |36) that grows exponentially as N —> oo . 

The coefficient c„„ and hence y m , can be expressed in terms of the root A,„. Indeed, from 
Eq.® we have 


1 

Cm 


n ( Am A ' m> )- a-a...~ p ' mUm) ■ 


i= 1 


^ /t /Im 


In particular, from the previous expression it is obvious that c m is real and positive, as it should 
be on account of Eq. ( [34] ). Another useful expression for the coefficient is obtained by noting 
that, by Eq. ( |30| >, 


m -1 

q m {A) = (A - 1 )p m (A) = d m+1 - 2 A' n + 1 =(A- \)(A - T m ) (d - d< m) ) , 

i=i 


and therefore 


1 


1 


lim 


q,„u ) q'mUm) (m + l)d™ - 2mA”; 


Cm Afu 1 A—>A m A A m A m 1 A m 1 

We thus arrive at the following remarkably simple expression for the coefficient y m in Eq. 


(37) 


A- 


ym 


(in + 1 )A m - 2m m — 1 

■ = m + 1 — 


C m 'l m 1 ^-m 1 

In particular, from Eqs. ( [3T| and ( [38] ) we obtain the following upper and lower bounds on A m : 

2m 


(38) 


m + 1 
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(39) 














which imply that 


lim A m = 2. 

m—*oo 

In fact, it can be easily shown that the sequence {/l m }“ =2 is monotonically increasing. Indeed, 
Pm +1 (y^m) Pm+lC^-m) ~ Pm(4 m ) ~~ —) ^ 0 ? 


so that A m < A m+ \ since A m+ \ is the largest root of the monic polynomial p m+ \. 

On the other hand, since A — 1 is not a root of p m , the characteristic equation p,„(A) = 0 can 
be written as 

lffl _ 1 

A m ~- -— = 0, (40) 


or equivalently 


Hence 


A- 1 


A = 2-. 

A m 


lim m(2 - A m ) = lim — = 0, 

m —>oo m —>cx3 si'" 


(41) 


and Eq. ([38]) implies that 


y m — 


A n/ f/l(2 A/n) 
A m - 1 


2 . 


Thus, for m » 1 we have 


2/V,min 


/ m \ N 
! 2 ’ 


m » 1 , 


where the RHS is the minimum average degeneracy of an su(p|^) model with p + q = m 
(cf. Eq. ([17])). 

It can be numerically checked that A m and (to a lesser extent) y m rapidly converge to 2 as m 
tends to oo (see, e.g., Table[l]below). To confirm this observation analytically, it suffices to note 
that from Eq. ([39]) we have 


1 1 

ira < 9 m 


1 + 


2 m ’ 


and therefore, by Eq. (|4T|), 


2 — — < A m < 2 . 

2>n 


The upper bound on A„, implies that y,„ < 2, while from the lower bound we easily obtain 

for m y 4 . 


m— 1 _ (m - l)e _ me 

y m > m +1-=2->2- 

1 - 2~"'e 2 m - e 2 m 


Remark 3. Equation ( |37] ) for the coefficient c m can in fact be extended with exactly the same 
proof to the remaining coefficients c- n) (i = 1 ), i.e.. 


(m) 

C. = 


a: 


(m) 


(drr>»+i)4 


(m) 
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Inserting the latter equation into Eq. ( [32] ) provides an alternative simple derivation of Binet’s 
formula for the w-nacci numbers mm 



Example 2. For m = 2 the roots of the characteristic polynomial /yl/f) are 


dt = ^(l-V5), A 2 = |(1 + V5), 


so that Eq. (|33)> reduces to the famous Binet formula for the Fibonacci numbers 




Thus the minimum average degeneracy |27} of a Yangian-invariant su(2) spin model is given by 



(42) 


in agreement with Eqs. ( |36| and ( |38[ i. 

In the case m = 3, the characteristic polynomial has two complex conjugate roots of modulus 
less than one and a real root 



Using Eq. ([38]) we easily obtain 



so that for large N the su(3) minimum average degeneracy behaves as 
d N ,min - 1.61702 • 1.63107^, N » 1. 


The exact formulas for /I 4 and y 4 are too unwieldy to display, and for m ^ 5 the characteristic 
equation defining A m cannot be solved in radicals. However, the numerical evaluation of A m and 
y m is straightforward; see Table[T]for a list of their values in the range m = 2 ,..., 10. 

4. Average degeneracy of translationally invariant models 

Several important examples of Yangian-invariant spin models possess the additional property 
of being translationally invariant. This is notably the case for the su(m\n) Haldane-Shastry spin 
chain (|3[-(|4j, whose energy function is of the form ( flO] ) with ;-;,v(./) = j(N - j). The dispersion 
relation of this chain obviously satisfies the identity 


enU) = e n (N - 
16 


(43) 








m 

A m 

Tm 

2 

1.61803 

1.38197 

3 

1.83929 

1.61702 

4 

1.92756 

1.76571 

5 

1.96595 

1.85899 

6 

1.98358 

1.91654 

7 

1.99196 

1.95139 

8 

1.99603 

1.97211 

9 

1.99803 

1.98420 

10 

1.99902 

1.99116 


Table 1: Values of A m and y m in the asymptotic formula J36} for m = 2,.... 10. 


characteristic of translationally invariant models (on a lattice with unit spacing between consec¬ 
utive sites) with a linear energy function. Indeed, for these models ey(./) is naturally interpreted 
as the energy of a quasi-particle with momentum 2nj/N. Since the momentum is defined up to 
multiples of 2n due to the translation invariance, the identity ( |43| > simply expresses the fact that 
the energy does not depend on the sign of the momentum. When the energy function En(8) is 
linear (cf. Eq. ([TO])) and its dispersion relation satisfies Eq. ( |43| ), the spectrum can be expressed 
as 

L(JV-1)/2J 

E N (8) = ^ s N (j)(6j + 5 N -j) + (1 - 7t(N))s n (N/ 2)6 N / 2 , 8 e A N (m\n), 

j= i 

where n(N) is the parity of N and |_-J denotes the integer part. Thus in this case E\ depends on 8 
only through the combination 8 + <i’ K , where <i’ K is the reversed motif with components ^ = S N .j. 

J J 

Conversely, if the energy function is linear and verifies 

E n (6) = f N (8 + 8 r ) (44) 


for some function fy, then the corresponding dispersion relation sn satisfies ( |43| . Indeed, it 
suffices to impose ( |44| ) for motifs of the form (0,..., 0,1,0,..., 0), which are allowed when 
m > 0, or of the form ( 1 ,..., 1 , 0 , 1 ,..., 1 ) when n > 0. Motivated by these observations, we 
shall say that a general Yangian-invariant spin model is translationally invariant whenever its 
energy function (not necessarily linear) satisfies the identity ( f44| >. The purpose of this section is 
to study the minimum average degeneracy of such models. 

By Eq. ( |44) . the minimum average degeneracy of a translationally (and Yangian) invariant 
spin model is given by 


j(s) 

Wmin 


(m + n) N 
v^\m\n) 


(45) 


where v N (m\n) is the cardinal of the set 


A§\m\n) = {5 + <S R | 8 e A N (m\n)}. 


(46) 


Given a motif 8 e Ay(/;i|n), we define its associated half-motif as the vector 


8 = 


(d>i + 8 n -\,82 + Sn-o, ■ ■ ■• 8 ( n - i )/2 + 8( jv+d/2), 

( 5 i + 8/v-i, 8 2 + 8n- 2, ■■ ■• 8(n-2)/2 + 8(n+2)/2, 8n/ 2), 
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N odd 
N even. 







Clearly, each vector 6 + 6 R is uniquely determined by its corresponding half-motif 6, so that 
in order to evaluate v^(m\n) it suffices to count the number of distinct half-motifs, ft is also 
important to observe that all the components of 6 can take the values 0,1,2 except 6 ^/2 = tf,v /2 
(when N is even), which can only take the values 0 or 1. 

In the genuinely supersymmetric case su(ra|n) with mu 4 0, it is straightforward to com¬ 
pute the number of distinct half-motifs, since in this case the motifs 6 are arbitrary elements of 
{0, 1 ) v ~ i . Thus, when N is odd the half-motif 6 is an arbitrary element of {0, 1,2} ,A l|/2 , and 
therefore 

v^(in\n) = 3 (A,_1 k 2 • mn ± q, jy odd. ( 47 ) 

On the other hand, when N is even we have 


Vjv (wi|n) = 2 • 3 {N 2 ; mn 4 0, N even. 


(48) 


since in this case the last component of 6 can only take the values 0 or 1. Hence the minimum 
average degeneracy of a translationally (and Yangian) invariant spin model in the genuinely 
supersymmetric case is given by 


j(s) 

^ At min 


V3 


m + n 


3 I m + n 

2 WT 


N odd 


N even. 


(49) 


(s') 

From the previous discussion it is obvious that u Nmm exactly coincides with the average degen¬ 
eracy d N when there is no accidental degeneracy, namely when the function f N in Eq. < |44[ > is 
injective. 


Example 3. Consider the su(m|n) elliptic model ( [281 in the special case m = n = 1, which is 
obviously invariant under integer translations. As shown in Ref. |36l , any spin chain of the form 


ff = 2*(i-/)(l- s S ll) )- 

i<j 


(50) 


where h is an even /V-periodic function, is isospectral to an su( 111) Yangian-invariant spin model 
with linear energy function and dispersion relation 


AM 

£/v(j) = 2(l-cos(2^7/iV))/t(/). 

/= l 


(51) 


Since a' v (./) = s(N - j), the dispersion relation © is also translationally invariant in the sense 
of the definition ( |44l >. In particular, for the su( 111) elliptic model ( |28| the dispersion relation can 
be computed in closed form using the techniques in Ref. Il36ll . with the result 


BnU) = 


N 

2 


( E\ N 2 1 

1 / \ 2 1 

l 1 KI + 4K 2 ^ 1 16 K’ 2 

P(j/N) - \£(j/N) -2rhjjJ 


(and sn( 0) = 0). In the previous formula K’ = V1 - k 2 J, 


E(k) 


= J' yjl - k 2 sin 2 ip d lip 

° 18 


j = l,... ,N — l (52) 












is the complete elliptic integral of the second kind with modulus k, p and £ denote the Weier- 
strass elliptic functions with period lattice generated by 1 and iK/(NK'), //] = £(1/2) and 
771 = £(1/2; l/2,i/VA''/(2A')). Using Eqs. ([TO} and ( [52] ) we have computed the spectrum of 
the su( 111) elliptic chain ([28} for N ^ 25 and several values of the modulus k. In all cases, the 
average degeneracy dw coincides with the minimum average degeneracy ([49} of a translationally 
(and Yangian) invariant spin model (cf. Fig. [4}. Thus the spectrum of the su( 111) elliptic chain 
does not exhibit any accidental degeneracy. 



N 

Figure 4: logarithmic plot of the average density of the su( 111) elliptic chain with k 2 = 1/2 (blue dots), compared to 
the minimum average degeneracy (49) of a translationally (and Yangian) invariant su( 111) spin model (continuous red 
lines, the top one for odd N and the bottom one for even N ). The dashed horizontal line represents the minimum average 
degeneracy 03 of a generic Yangian-invariant su(l| 1) model. 


In the non-supersymmetric (purely bosonic or fermionic) case, the computation of is 
considerably more involved, due to the fact that not all elements of {0, 1 )' v 1 are valid motifs. 
For this reason, we shall restrict ourselves in the rest of this section to the simplest and most 
important su(2) case. In the bosonic case (which, as we know, is equivalent to the fermionic 
one), an su(2) motif 6 e {0, 1 }' v_l cannot contain two consecutive l’s. Hence, for odd N the 
corresponding half-motif 6 e {0,1, 2\ iN ~ Y ^ 2 satisfies the following two constraints: 

i) No 2 can be preceded or followed by a 1 or a 2. 

ii) The last component of 6 cannot be equal to 2. 

It is also easy to see that any element of {0,1, 2\ (N ~ y ^ 2 satisfying these conditions is a valid half¬ 
motif 6 when N is odd. On the other hand, when N is even the last component of a half-motif 
is the middle component of the corresponding motif, and therefore condition ii) above should be 
replaced by the more stringent one 

ii') The half-motif 6 can only end in 0 or 01. 

As before, it can be shown that any element of {0,1, 2} N ^ 2 satisfying conditions i) and ii') is a valid 
half-motif. Thus, in order to compute the minimum average degeneracy for odd (respectively 
even) N we need only count the number of elements of {0,(respectively {0,1,2}^ 2 ) 
satisfying conditions i) and ii) (respectively i) and ii')). 
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Consider, to begin with, the case of odd N - 2r + 1, and denote by p, the number of vectors 
6 e {0, 1,2) r satisfying the constraints i) and ii). Clearly, any such vector 6 can only end in xl, 
jrO and 20, where x is either 0 or 1. In the first two cases, the first r — 1 components of 6 still 
satisfy conditions i) and ii), so that the number of valid half-motifs of length r ending in xl or 
xO is equal to 2p r _i. On the other hand, if 6 ends in 20 its first r — 2 components are a vector 
belonging to the set of half-motifs of length r — 2 not ending in 1 , whose cardinal is obviously 
H r _2 - Hr- 3 . We thus arrive at the recursion relation 

Hr = 2flr-\ + Hr- 2 ~ Hr-3 ■ (53a) 

It is straightforward to check that h\ = 2, hi =5 and H 3 = 11, which, together with ( |53a| ), 
completely determines Hr for all r. In fact, it is immediate to verify that the latter initial conditions 
are equivalent to the simpler ones 


H-2 = H-i - 0, Ho- 1- (53b) 

Consider next the case of even N = 2r, and let fi r denote the number of vectors 6 e {0, 1,2) r 
satisfying the constraints i) and ii')- The reader should easily convince himself that the valid 
motifs are precisely those from the IV = 2r + 1 case not ending in 11, so that 

Hr = Hr - Hr—2 ■ (54) 

Since Hr (and hence Hr - 2 ) is a solution of the linear recursion relation with constant coeffi¬ 
cients ( |53a| ), it follows that jl r satisfies the same recursion relation: 

Hr = 2/j,._i + Hr- 2 ~ Hr —3 ■ (55a) 

However, the initial conditions are now h\ - 2, /H = 4, jj\ - 9, or equivalently 

H-2 = H-\ = Ho = 1 - (55b) 


Thus, the minimum average degeneracy of an su(2) translationally (and Yangian) invariant spin 
model is given by 


j(s) 

^ At min 


v^( 2) 


(56) 


where 


v^(2) = 


H(N- 1)/2 , 
HN /2 , 


N odd 
N even. 


(57) 


The characteristic polynomial p(A) = A 3 - 2 A 2 -4+1 of the recursion relation d53a)-(|55al 
has three real roots given by 



V7 cos(i/j + 2kn/3)^ , 


k = 0,1,2, 


with 


= \ arctan(3 V3). 
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As explained in Section 3.2 the sequences Hr and fi r (r = 1,2,...) can be expressed as 

2 2 

Hr — > Hr ~ Ok'\ » 

k= 0 t=0 

where the coefficients a * and d /( are respectively determined by the initial conditions ( |53b| ) and 
( |55b[ ). (In fact, from Eq. ( [54] ) it follows that a* = a*(l - d7 2 ).j Since the root of p(A) with largest 
absolute value is do - 2.24698 (while [/I12I < 1), the asymptotic behavior of Hr and fi r as r —» 00 
is given by 

Hr — a<) 4q , Hr — ao Aq . (58) 

r —>00 r —>00 

The coefficients ao and do can be readily computed in closed form, with the result 
4(1 + V7 cos^) 2 ^ - 


ao 


21(1 + 2 cos 2<p) 


4(1 + V7 cos w) - 9 

0.97869, do = - -—-=* 0.78485 . (59) 

21(1 +2 cos lip) 


From Eqs. ( ]57| > and ( |58] > we easily deduce the asymptotic behavior of the minimum average 
degeneracy dj^^ as N —> 00 in the su(2) case, namely 


V,min 


r“> 1—f 

,! -"W ' 


(60) 


where the coefficient is given by 


Vd0 


(s) 

Tyv = 


ao 

1 

a 0 ’ 


N odd 
N even. 


(61) 


As N —> 00, this minimum average degeneracy becomes considerably larger than the correspond¬ 
ing one for a generic Yangian-invariant su(2) model. Indeed, from Eqs. (|42[> and d6Q[) we obtain 


i (s) 

iV.min 


dN ,min N—>oo 


- a N 


(V5 - 1) Vdo 


=* a N x 1.07941 - N , 


where a N = 2y? L/(5 - V5) is close to 1 (1.10830 for odd N and 0.92197 for even N). 


5. Average degeneracy of spin chains of Haldane-Shastry type 

Ever since Haldane’s original paper 0, it has been noted that the spectrum of the HS chain 
exhibits a high degree of degeneracy, which is generally attributed to the Yangian symmetry of 
this model. In fact, since both the su(n?|n) HS and PF chains are Yangian-invariant spin models, 
their average degeneracy is bounded below by du.mm in Eqs. l | 1 A\ (in the genuinely supersymmet¬ 
ric case) or ( fTTj i (in the purely bosonic or fermionic case). In the case of the FI chain, although 
its symmetry under the Yangian has not been proved rigorously, since the spectrum is given by 
Eqs. 0 and 0 the same bounds apply. Roughly speaking, for the simplest su(2) case it can 
be checked that the average degeneracy of the HS, PF and FI chains is of the same order of 
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magnitude as the above bounds —or, in the case of the HS chain, their refinement (|55)-((37) for 
translationally invariant spin models— only for low values of N (say, N < 10). In this sec¬ 
tion we shall study the average degeneracy of these models for much higher values of N (up to 
N = 50), and show that, contrary to widespread belief, it cannot be solely explained by their 
Yangian invariance. 

As explained above, for simplicity’s sake we shall restrict ourselves to the su(2) (bosonic 
or fermionic) case. For the PF chain, the computation of the average degeneracy is completely 
straightforward, since its spectrum is a set of consecutive integers lf39l and the minimum and 
maximum energies can be found without difficulty. Indeed, consider (for instance) the bosonic 
case. The minimum energy is clearly zero, corresponding to the su(2) motif 6 = (0,.... 0). On 
the other hand, Eqs. © and m imply that the maximum energy is obtained from the motif 
(..., 0,1,0,1), with the result E Nt max = ( N 2 - tt(N))/4. Thus the number of distinct energy levels 
is given by 

N 2 - n(N) 

V N (2) = F At max + 1 = -^-- + 1 (PF chain), (62) 

and the average degeneracy is of course c/ v = 2 N /v N (2). In the case of the HS or FI chains, 
no simple way of expressing in closed form the number of distinct levels as a function of N is 
known. For both of these chains, the most efficient way of computing the spectrum for relatively 
high values of N is to evaluate the partition function lfl9l(4Ttt in a recursive way. More precisely, 
it can be shown that the partition function Z ;V of the (fermionic) FI chain is determined by the 
recursion relation 


Z N+l (q) = 2+ (1 - q N ^ N - l W N - l)(a+N ' 2) Z N -M\ q - e" 1 '^ , 

with the initial conditions Z_\{q) = 0, Z {) (q) = 1. The situation is a bit more involved for the 
HS chain, since its dispersion relation |6]i depends explicitly on N. In this case, the (fermionic) 
partition function is given by Z v (<y) = Q,\:(q\ N), where Q k (q',N) is obtained from the recursion 
relation 

Q k +i(q',N) = 2q k(N ~ k) Q k (q; N) + (1 - 

with the initial conditions Q-\(q\N) = 0, Qoiq'-N) = 1. Proceeding in this way it is feasible 
to compute the partition functions of the su(2) FI and HS chains for N ^ 50 in a standard 
desktop computer, and thus determine the number of distinct energies and the corresponding 
average degeneracy d^. In Figure [5] we present logarithmic plots of dn for all three chains of 
Haldane-Shastry type with N ranging from 10 to 50, compared to its lower bounds ( [27] ) (for 
the PF and FI chains) and (|56|)-(|57j) (for the HS chain). It is apparent from this figure that as N 
becomes moderately large (of the order of 25) d^ grows much faster than expected for a generic 
Yangian-invariant spin model. The ultimate reason for this behavior is the fact that these chains 
have a linear energy function whose dispersion relation is a low degree polynomial. Indeed, we 
shall show that this implies that the number of distinct levels has at most polynomial growth 
with N, compared to the typical exponential growth for a generic Yangian-invariant spin model; 
cf. Eqs. ( f]~6l >, (|26])-(f34]>. (|47|)-(|48|), (|57]>-(|58]>. As a matter of fact, we have already proved that the 
number of distinct levels of the su(2) PF chain is a second degree polynomial in N; see Eq. ( |62[ i. 
We shall next prove an analogous result for the su(2) HS and FI chains. 

In the case of the HS chain, since the dispersion relation g,v(./) is a polynomial in j with 
integer coefficients, and e,v(y) > 0 for 1 ^ / ^ <V - 1, the energies are clearly nonnegative 
integers. Moreover, from Eq. ([6] it follows that all the energies are even when N is odd. Since 
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Figure 5: Left: logarithmic plot of the average degeneracy of the Haldane-Shastry su(2) spin chain (blue circles), com¬ 
pared to the minimum average degeneracy of a translationally (and Yangian) invariant spin model (red circles). Right: 
analogous plot for the su(2) Polychronakos-Frahm (blue circles) and Frahm-Inozemtsev (blue triangles for a = 3, blue 
squares for a irrational) spin chains, versus the minimum average degeneracy of a generic Yangian-invariant su(2) spin 
model (red circles) and of an su(2)-invariant spin chain of the form (3) (gray circles). 


the minimum energy in the bosonic case is zero, the number of distinct energies satisfies the 
inequality 


v n (2) < 


E iV,max ^ j 

1 + n(N) + ’ 


where E ^,max is the maximum energy of the bosonic chain. By Eqs. (|22[i-(|23[i. this maximum 
energy is given by 


-'Y,max 




: (F) 

Y,min 


7=1 


1 

6 


N(N 2 - 1 )-E 


(F) 

Y,min ’ 


where 

<L = ^(iV 2 -4 + 37r(iV)). 

is the minimum energy of the fermionic chain computed in Ref. OTj . We thus obtain 


vw(2) < • 


N , 

— (V 2 +2)+l, 
N , 

i24^- 1) + 1 ’ 


N even 

N odd 


(HS chain). 


(63) 


In fact, as can be seen from Fig. [6] these bounds provide a good approximation to v,y(2) for 
large N. 

Finally, in the case of the FI chain, the behavior of v N (2) is quite different when a is rational 
or irrational. Indeed, if a = a/b is an irreducible rational number then bE N (6) is a nonnegative 
integer. Thus the number of distinct energies is bounded above by bE^, max + 1, where £)v,max is 
again the maximum energy of the bosonic chain. Proceeding as before and using the value of 
Inin com P u ted in Ref. ED we readily obtain 


vn( 2) < • 


N 

12 

1 


V 12 


(2 bN 2 + 3 (a - b)N -2b) + \, 
(N 2 - 1)(2 bN + 3 (a - b)) + 1, 


N even 


(FI chain, a e Q). (64) 
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Figure 6: logarithmic plot of the number of distinct energy levels of the su(2) Haldane-Shastry chain (red circles for even 
N, blue circles for odd AO, compared to the upper bounds {63} (solid lines). Inset: same plot for the range 30 ^ N ^ 50, 
in ordinary (non-logarithmic) scale. 


On the other hand, when a is irrational we can write 

N-\ N -1 

E n ( 6) = jdj + 2 j(J - 1)6j = aE ( °\6) + E%\6), 

7=1 7=1 

where the functions E { ^(6) (i = 0,1) take nonnegative integer values and (in the bosonic case) 
actually vanish when 6 is the zero motif. Hence 


vtv( 2 )<(£® +l)(£jvl x + l), 


a) 


"'A^max 


Af,max 


where £^ max is the maximum value of E ( !j(6) over all valid su(2) bosonic motifs 6. Although 
Ovm-K can be computed without difficulty, for our purposes it suffices to note that 


Tv, max 


IX = ? (N-1), <L<Zi(;-i) = f (iV-l)(JV-2), 

7=1 7=1 

from which it easily follows that 

N 5 


vn( 2) ^ 


(FI chain, a i 1 


(65) 


Although the previous discussion has been restricted to the su(2) Haldane-Shastry type 
chains, its main conclusions extend to the whole class of su(m|«) Yangian-invariant spin models 
studied in Refs. ED ED, whose energy function is of the form ( fT0| ) with a dispersion relation 
£ v(./) polynomial in j and N. Indeed, suppose that 


svO') = X a rs N''f , 


( 66 ) 


( r,s)el 


where I is a finite subset of No x No, No being the set of nonnegative integers, and a rs is inde¬ 
pendent of N and non-zero for all (r, s) 6 I. The energy function is then given by 


TV—1 


E n (6)= J] »V,.Y' X./"<>/ = X a r S N r E^(8), 


(r,s)el 


j= 1 
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( r,s)el 





where each function E^\S) takes nonnegative integer values and is bounded above by the power 
sum 

1WS+1 

S^= Yj' = -- + 0(N S ). (67) 

s + 1 

7=1 

We thus obtain the polynomial bound 

v N (>n\n) < Y] (S« + 1) = Y] + l) rW , (68) 

(r,s)el s 

where r(s) is the number of indices r such that (r, s ) e /, and the last product ranges over all ,v 
such that r(s) > 0. Note that, on account of Eq. ( |67] , the bound ( [68] ) satisfies 

Y\ (S“ + l) r(s) = o(N Zs(j+1,rW ) . 

S 

The bound in Eq. ( |68] > can be considerably strengthened when all the coefficients a rs in the 
dispersion relation ( |68] > are rational numbers, as is the case with the PF and HS chains (or the 
FI chain, when a e Q). Indeed, assume that a rs = a rs /b, where a rs e Z and b e N. In this case 
bE^iS) takes only integer values, so that 

Vy(m|u) ^ b^Eft ' max — E.Y^iin) + 1, 

where E’Y.min an d E Nmax denote respectively the minimum and maximum energies. From the 
inequality 

M/)| ^ 2 \ a ^ Nr ] S < N " Z |ar "! = ANk . 

(r,s)el ( r,s)el 

where k = max(r + s ) is the total degree of t'N(j) as a polynomial in ( j,N ), we immediately 

(r,s)el 

obtain 

En, max - E n% min < 2 max |^(d)| ^ 2N max \s N (j)\ ^ 2 AN k+l , 

6eA N (m\n)' 1 1 1 

and hence 

v N (m\n) < (2bA + 1 )N k+1 . (69) 

Note that this equation is in agreement with the asymptotic behavior of the bounds ([62])—([64]) for 
the PF, HS and FI chains (with a e Q in the latter case). 

6. Anyons and Fibonacci numbers 

The elementary excitations of certain quantum many-body systems can be described in terms 
of anyons , i.e., effective quasi-particles obeying neither the Bose-Einstein nor the Fermi-Dirac 
statistics ESEl. The original definition of anyons is intrinsically two-dimensional, since it 
requires that the wave function change by a suitable phase under exchange of two identical 
particles, which is only consistent in two dimensions. An alternative definition of anyons in 
terms of a generalized exclusion principle was proposed by Haldane lH4l . According to this 
definition, in a system of (identical, non-interacting) anyons the number of available one-particle 
states (or orbitals) decreases by gAk when the number of particles increases by A k. The constant 
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g, which must be a rational number, is called Haldane’s exclusion parameter. Clearly, bosons and 
fermions have exclusion parameter g = 0 and g = 1, respectively. The advantage of Haldane’s 
definition is that, in contrast with the original one, it is also valid in one dimension. For instance, 
it is well known that when the coupling constant is a nonnegative integer a the spectrum of the 
/V-particle Calogero model ll29l is that of a system of N free anyons with exclusion parameter 
g - a in a harmonic oscillator potential (see, e.g., the review paper |48|). 

Another important example of one-dimensional anyons is furnished by the su(2) Haldane- 
Shastry chain, whose “spinon” excitations (in the antiferromagnetic case) are known to be anyons 
with g = 1/2 lfl4l . It should be clear that this property is in fact shared by all su(2) (fermionic or 
bosonic) Yangian-invariant spin models, regardless of their energy function. To see this, consider 
an su(2) (say, bosonic) motif 6 and its associated border strip k = (k\ ,..., k,). where k\ + ■ ■ ■ +k r = 
N is the number of spins. In the su(2) case, it is a simple matter to compute the dimension 
of the corresponding T(gl(2|0))-module 14(2|0) = Vf. Indeed, a (bosonic) su(2) border strip 
k cannot contain rows with more than two boxes. Moreover, any tableau corresponding to k 
must necessarily contain the sequence 01 (from left to right) in any row with two boxes. The 
dimension of 14 is thus given by the number of ways of filling the remaining (“empty”) boxes in 
each column of the border strip k (from top to bottom) with sequences of the form 0 ■ ■ ■ 01 ■ • • 1 
(where the number of 0’s or l’s can of course be zero). Thus the r columns of k play the role of 
available orbitals, and its N - 2r + 2 empty boxes that of particles with two internal degrees of 
freedom (0 or 1). (These particles are obviously bosons, since each orbital can be occupied by 
more than one particle with the same quantum numbers.) Thus in this case we have k = N-2r+2 
and A r = -Ak/2, so that Haldane’s definition is satisfied with g = 1/2. This shows that any su(2) 
Yangian-invariant model is indeed equivalent to a system of anyons with two internal degrees of 
freedom and exclusion parameter g = 1/2 ( semions , in Haldane’s terminology). 

Haldane’s definition of anyons can be equivalently formulated in terms of the total number 
of states in a system with k anyons and N one-particle states (or orbitals). Indeed, for anyons 
with exclusion parameter g this number is given by 



(70) 


When g is a natural number, a very simple realization of this statistics can be obtained by consid¬ 
ering a system of k particles and N orbitals satisfying the following generalized Pauli exclusion 
principle: two occupied orbitals must be separated by at least g - 1 empty ones. This requirement 
is clearly reminiscent of the rule satisfied by the bosonic su(2) motifs. More precisely, we can 
interpret such a motif 6 e A,v(2|0) as encoding a state of a quantum system with N - 1 orbitals 
(the components of 6), with a 1 indicating an occupied orbital and 0 an empty one. With this 
interpretation, the 1 ’s in the motif 6 are quasi-particles satisfying the generalized Pauli principle 
with g = 2, and thus can be regarded as anyons with exclusion parameter g = 2. This point 
of view is particularly natural when the energy function is linear (cf. Eq. ©), as is the case 
with spin chains of Haldane-Shastry type (cf. Eqs. <0-Q). Indeed, for such Yangian-invariant 
models the energy function can be written as 


r -1 


E N (fi) = YjSNiKi), 


(71) 


where Kj is the position of the z'-th 1 in 6. It should be stressed that this interpretation is not 
equivalent to the one discussed above in terms of spinons, for which g = 1/2. Indeed, in the 
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latter interpretation the l’s in the motif (or, more accurately, the sequences (0,1)) simply separate 
the orbitals, while the remaining 0’s correspond to the particles. In particular, a fundamental 
difference between both interpretations is that in the new one the degeneracy associated to each 
motif (that is, the dimension of the corresponding invariant module VjJ is not taken into account. 
For this reason, the new interpretation defines a different type of system, whose Hilbert space 
dimension equals the number v^(2) of distinct su(2) motifs. Since v,v(2) is given by the Fibonacci 
number F n+ i (cf. Eq. ([26])), from Eq. ( [70] ) with N — N - 1 and g - 2 we obtain the remarkable 
identity 



where the upper limit is determined by the condition N - k ^ k. In other words, the (N + l)-th 
Fibonacci number is obtained by diagonally summing the binomial coefficients arranged in a 
Pascal triangle lf34l . This classical number-theoretic result can thus be interpreted as expressing 
the dimension of the Hilbert space of a system of g — 2 anyons with N - 1 orbitals in terms of 
the corresponding dimensions of its k-particle subspaces. 

Motivated by the previous discussion, we shall next consider the bosonic su (m) motifs with 
m > 2 as representing the states of a system of quasi-particles (namely, the l’s in the motif) 
with N - 1 orbitals (the motif’s components), subject to the exclusion principle that there can 
be no more than m - 1 consecutive occupied orbitals. As before, this rule obviously does not 
take into account the degeneracy associated to each motif, so that the system under consideration 
is not an su(m) Yangian-invariant spin model. It can instead be regarded as an effective model, 
whose spectrum ( |7T] > reproduces the set of distinct energy levels of a generic (bosonic) su(m) 
Yangian-invariant spin model with linear energy function and dispersion relation s^ij ). 

We shall next compute the statistical weight w* for the new system just introduced. To this 
end, let 6 e Aai(w| 0) be a bosonic su(m) motif with jo single l’s (i.e., (1,0) sequences), j\ double 
l’s (i.e, (1,1,0) sequences), ..., and j m -2 sequences (1,..., 1,0), with 

m -1 

jo +2ji + ■■■ + (m - 1 )j m -2 = k . (73) 


Removing the zero in each sequence (1,0), (1,0) in each sequence (1,1,0), ..., (1,..., 1,0) in 

m-2 

each sequence (1,..., 1,0), except for the zero (if any) following the last 1 in the motif, we are 

m -1 

left with a sequence 6 e {0,1 } N ~ k containing jo + ■■■ + j m -2 l’s. Conversely, it is clear that 
from any such vector 6 e {0, l}^ - * one can construct a bosonic su (m) motif with jo single l’s, j) 
double l’s, ..., and /,„_2 sequences (1,..., 1,0) by adding a zero after jo l’s, the sequence (1,0) 


m -1 


after j\ of the remaining l’s, etc. Thus, the number w(j 0 ,..., jm- 2 ) of motifs with jo single l’s, 
/j double l’s, etc., is given by the number of ways of choosing the positions of the jo l’s starting 
a sequence (1,0), the j\ l’s starting a sequence (1,1,0), etc., among the integers 1 — In 

other words, 


w(j 0 , ■ ■ •, jm-2 ) 


(N - k)\ 

m-2 m-2 


w-k-z jj )j n u\ 

i =0 1=0 
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The statistical weight w* is the sum of all w(jo ,..., y m _ 2 ) whose arguments (jo,... ,j m _ 2 ) sat¬ 
isfy f73} , namely 

V ( N-k)l 

w k = 2j -—-ITT— ■ (74) 


, m—2 x nt 

JO,~,U-2>0 (N-k- V /■■)! n 7«! 

7o+2; 1 +-+(m-l); m _ 2 =* V iff, J, > Jl 


m—2 

n 

i=0 


Note that for m = 2 this equation obviously reduces to Eq. CD with g = 2 and N — N - 1. 
Eliminating jo with the help of Eq. ( |73} and explicitly enforcing the condition 2”i 0 2 ji ^ N - k 
arising from the first factor in the denominator we obtain the equivalent expression 


w k = 


E 


(N - k)\ 


2y'i +3 ji +• • • +(m- 1 )j m - 2 ^ k 
ji + 2 / 2 +—+ (m- 2 )jm -2 ^ 2k-N 


m -2 , . m—z v m-z 

(N - 2k + y EO'+D*V n 7«! 

i= 1 7 v -- 1 7 -- 1 


m—2 

E< 

i=l 


m-2 

n 

i=l 


(75) 


For instance, for m = 3 we have 

L*/2J 




Z 


(N - k)\ 


j=max(2k-N,0) 


(N — 2k + j)\(k - 2j)\j\ ’ 


m = 3. 


The total dimension of the system’s Hilbert space, which is the same as the number v N (m) of 
bosonic su(m) motifs for a Yangian-invariant spin model with N spins, is obtained by summing 
the statistical weight w* over all possible values of k. If N - 1 = mi + s with 0 y s ^ m - 1, the 
maximum value of k is clearly 


, ,, m — 1 , m - 1 m — 1 — s 

r(m — 1) + s = - (N — 1 — s) + s = - N - 


m - 1 


■ N 


since O^m-l-s^m-l. From Eqs. ([75} and ([26} we finally obtain the identity 


1 N+m-l 


\_{m-\)N I m\ 

z 

k =0 


z 


(N - k)\ 


Jl .-Jm-2^0 
2y i +3 ji +—+(m-1 )j m -2 ^ & 

7 1 +2^2+• • • +(m-2)j m -2 ^ 2 *-W 


m-2 


m-2 


m—2 


(N-2k+ y y'/)!(*- 2 0+1)7'.-)! n 7'i! 

i'= 1 /= 1 i= 1 


(76) 


An alternative version of the latter equation follows by using Eq. ( [74] ) for Wk and eliminating k 
using Eq. CD- namely 


F (/») 

1 N+m -1 


z 

70vJm-2^0 

2 y'o+3 7 i H — \-mj m -2 


m-2 

A- y (i + 1)7/)! 

/=0 

. m—2 . m-2 

Ov- E 0 + 2)7,-)! n 7i! 

v i=0 1=0 


(77) 


This is essentially the classical combinatorial expression for the generalized Fibonacci numbers 
proved in lf34l (cf. Eq. (28") in the latter reference). 

It is apparent from Eqs. ( [70} (with N = N- 1) and ( |75} that for m > 2 the quasi-particles asso¬ 
ciated to the (bosonic) su(m) motifs of a Yangian-invariant spin model are not anyons. Note, how¬ 
ever, that there is a generalization of Haldane’s definition ( [70} , due to Murthy and Shankar 149] , 
to systems with an infinite-dimensional one-particle Hilbert space. The latter generalization is 
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essentially based on the behavior of the statistical weight (J70]i as the number of orbitals N tends 
to infinity for any fixed k, namely 




(78) 


from which it follows that 


1 N 

2~ s = !klw 


v (k\w k (g) \ 
-1)1 N k r 


(79) 


This observation suggests looking at the asymptotic behavior of the density of states CD as 
N — N - 1 —* oo with k fixed. To this end, note first of all that j, ^ k/(i + 1) on account of the 
restriction X™“j 2 (; + 1 )j, ^ k, and that 


(N - k + 1)! 

(N — 2k + 1 + 1)! jv~*oo 


N k -' 


+ 0(N k -'~ l ). 


Thus the leading terms in the sum ( |75) > are obtained for j\ = ■ ■ ■ = j m -i = 0 and j\ = 1, 
j 2 = • • • = jm -2 = 0, or more precisely 


Wk 


(N - k + 1)! 
(N -2k+\)\k\ 


(N-k+l)\ 

+ (N -2k + 2)\(k-2)\ 


+ O (N k ~ 2 ) 


lN-k+l\ N k ~' 

l k ) + (k-2)\ 


+ 0(N k - 2 ). 


Using Eqs (f70|-([78| with g - 2 to expand the first term on the RHS we immediately obtain 


N k 


Wk= ld 


1 N k - X 

2 (k-2)\ 


+ O (N k - 2 ). 


Thus as TV — > oo with k fixed the density of states ( |75[ i satisfies Eq. ( [78] > with g — 1. Note that, as 
already observed by Murthy and Shankar, if the system’s maximum and minimum energies are 
independent of N (or, more generally, the spectrum is bounded by A/'-independent constants) the 
limit N —> oo with k fixed can be naturally interpreted as the continuum limit. We conclude that 
in the continuum limit the (bosonic) su(m) motifs with m > 2 behave as a system of fermions, 
according to Murthy and Shankar’s definition. 


7. Conclusions 

In this paper we consider a wide class of finite-dimensional Yangian-invariant quantum mod¬ 
els —which we have called Yangian-invariant spin models— closely related to representations 
of the Yangian F(gl(m|n)) labeled by border strips. This class includes several widely studied 
systems, such as spin chains of Haldane-Shastry type and certain families of one-dimensional 
vertex models. A common feature of all Yangian-invariant spin models is the relatively high 
degeneracy of their spectrum, essentially due to their symmetry under the Yangian. As a quan¬ 
titative measure of the degree of degeneracy of these models, we compute in closed form the 
intrinsic average degeneracy of any such model stemming from its Yangian invariance. It turns 
out that, in the non-supersymmetric case, this minimum average degeneracy is expressed in terms 
of generalized Fibonacci numbers. Using several fundamental properties of these numbers, we 
derive the asymptotic behavior of the minimum average degeneracy as the number of spins tends 
to infinity. Our results provide a stringent test for a quantum system to be a Yangian-invariant 
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spin model. As an example, we apply this test to the su(»;|n) supersymmetric version of Inozemt¬ 
sev chain, which has been studied in the context of the gauge-string duality as a candidate for 
reproducing the spectrum of the dilation operator at several loops; see, e.g., 123. We have also 
refined the previous results on the minimum average degeneracy for models that are in addition 
translationally invariant, both in the su(»i|n)-supersymmetric and in the simplest su(2) cases. 

The second main problem that we address in this paper is the analysis of the precise role 
played by the Yangian symmetry in the high degeneracy of the spectra of spin chains of Haldane- 
Shastry type. We show that these chains are much more degenerate than generic Yangian- 
invariant spin models, since the number of distinct levels grows polynomially for the former 
systems and exponentially for the latter. This accidental degeneracy is essentially due to the 
polynomial character of the dispersion relation, and is therefore shared by the family of one¬ 
dimensional vertex models studied in Refs. mm. Moreover, if the coefficients in the disper¬ 
sion relation are rational numbers (as is the case for the HS and PF chains), we show that the 
accidental degeneracy is even higher, as illustrated by the FI chain with a rational vs. irrational 
parameter. Finally, we consider an effective model of quasi-particles describing the distinct en¬ 
ergy levels of a non-supersymmetric Yangian-invariant spin model. We show that in the su(2) 
case these quasi-particles are anyons with exclusion parameter g = 2, while in the su(m) case 
with m > 2 they become fermions in the continuum limit. We also provide a simple interpreta¬ 
tion of several combinatorial expressions for ra-nacci numbers in terms of the statistical weights 
of the effective model of su(m) quasi-particles. 
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